Resonant normal form for even periodic FPU 

chains 



Andreas Henrici* Thomas Kappeler^ 
February 1, 2008 

Abstract 

In this paper we investigate periodic FPU chains with an even number 
of particles. We show that near the equilibrium point, any such chain 
admits a resonant Birkhoff normal form of order four which is completely 
integrable - an important fact which helps explain the numerical experi- 
ments of Fermi, Pasta, and Ulam. We analyze the moment map of the 
integrable approximation of an even FPU chain. Unlike in the case of 
odd FPU chains these integrable systems (generically) exhibit hyperbolic 
dynamics. As an application we prove that any FPU chain with Dirichlet 
boundary conditions admits a Birkhoff normal form up to order four and 
show that a KAM theorem applies 

1 Introduction 

In this paper we consider FPU chains with an even number TV of particles of 
equal mass, normalized to be one. Such chains have been introduced by Fermi, 
Pasta, and Ulam [5] , as models to test numerically the principle of thermaliza- 
tion as N grows larger and larger. A FPU chain consists of a string of particles 
moving on the line or the circle interacting only with their nearest neighbors 
through nonlinear springs. Its Hamiltonian is given by 

N N 
n— 1 n— 1 

where V : M. — > R is a smooth potential. The corresponding Hamiltonian 
equations read (1 < n < N) 

in = d Pn Hv = p n , 

pn = -d qn Hv = -V'(q n - q n +i) + V'(q n -i ~ q n )- 
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Here q n denotes the displacement of the n'th particle from its equilibrium 
position and p n is its momentum. If not stated otherwise, we assume periodic 
boundary conditions 

(q i+N ,p i+N ) = (qi,Pi) Vi S {0, 1}. (2) 

Without loss of generality, the potential V : K — > K is assumed to have a 
Taylor expansion at of the form 

V(x) = k fa - |* 3 + ^ + . . .) , (3) 

where k is the (linear) spring constant normalized to be 1 and a, (3 € K are 
parameters measuring the strength of the nonlinear interaction. The minus sign 
in front of the parameter a in the expansion ([3]) turns out to be convenient 
for later computations. Substituting the expression §5§ for V into (fl]), the 
corresponding expansion of Hy is given by 

I N I N a N r N 

Hv= 2^ P " + 2 H^+rtn)^ H(*H-rtn) 3 +^ ^(<Z«frtn) 4 +- ■ ■ ■ (4) 

n— 1 n— 1 n— 1 n— 1 

For any FPU chain, the total momentum P = -4° Xm=i ^™ ^ s an integral of 
motion, and therefore the center of mass Q — -h ^] 1 g n evolves with constant 
velocity. Hence any FPU chain can be viewed as a family of Hamiltonian systems 
of 2N — 2 degrees of freedom, parametrized by the vector of initial conditions 
(Q,P) S K 2 with Hamiltonian independent of Q. In particular, for N = 2 
any FPU chain is integrable, and hence we will concentrate on the case N > 
3. Further note that for any vector (Q,P) € R 2 , the origin in R^" 2 is an 
equilibrium point of the corresponding system. The momentum of such an 
equilibrium point is given by the constant vector (pi, . . . ,Pn) = P (1, . . . ,1). 

The frequencies (a^)i<fc<jv_i of the linearization of an arbitrary FPU chain 
at (q,p) = (0, 0) can easily computed to be 

o n ■ k 7 *' 
uj h = 2 sin — . 

k N 

The corresponding resonance lattice is given by (see Appendix A of [7]) 

jz = (h, . . . , i N _!) e z^- 1 1 k sin 77 = °| 

and generated by the vectors l^- k \ 1 < k < N — 1, defined by — ek — £N-k, 
where e,-, 1 < i < N — 1, denotes the standard basis in l^ -1 . 



For any point (ir,y) S R introduce the variables / = (/fe) 



Kfc<jV-l, 



where 



h = \{xl + yl), (5) 
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as well as M = {M k )x<k<N-i, J = (Jk)i<k<N-i, and L = (Lfc)i<fc<iV-i- They 
arc defined on R 27V ~ 2 with values in R^ 1 and given by 

M k ^^(xkUN-k-XN-kyk) ; Jk = ]^{xkXN-k+VkVN-k) ; L k = -{l k - i N - k ) • (6) 

Note that for any 1 < k < N - 1, (M k , J k ,L k ) = (-M N _ k , JN-k, ~LN-k) and 
I h I N -k = Ml + Jl, or 



Ik + IN-k 



M 2 + J 2 + L 2 , 



(7) 



i.e. Mfe, J kl L k are the Hopf variables expressed in x k , y k , xjv-fc, yjv-fc- They 
describe the image of the Hopf map from the three-dimensional sphere of radius 
+ IN-k) centered at the orgin of R 4 . Further define the function H a ^ : 
I, given by 



pAf-l 



N-l 



N-l 



, , v-^ kn 1 v-^ . 9 8— a v-^ lit m7r , n 

Ha Ai) = 2J2 ™ ^ h + -fit E d Pl + -^r E sin n sm at 1 * 7 ™ (8) 

l<l,m<N-l 



fe=l 



fc=l 



where 



and let 



where 



and 



A:tt 



d+ = d+(a,8) :=a 2 + (/3-a 2 )sin 2 — 



R a , (J,M) := 



P -a 2 
AN 



(r(J,M) + Rk(J,M)) 



R(J,M) = 4 sm^(j fc J f _ fe -M fc M f 



i<fc<^ 



i?™ (J, M) 



7 2 - M 2 

4 4 





JV 
4 

otherwise. 



if f e N 



(9) 
(10) 
(11) 

(12) 



Note that for a, 8 £ K with /3 = a 2 - referred to as Toda case -, R a .j3 vanishes. 

In [7] we proved that even FPU chains admit a resonant normal form near 
the equilibrium point. 

Theorem 1.1. If N>A is even, there are canonical coordinates (x k ,y k )i<k<N-i 
so that the Hamiltonian of any FPU chain, when expressed in these coordinates, 
takes the form H# unc (I, J, M) + 0(\(x, y)\ 5 ) where 



Hi 



NP 2 



H a , p (I)-R a ^(J,M) 



(13) 



with a, 8 £ K as in (0) and where H a ^(I) and R a ^{J, M) are given by (0j and 
HO)) , respectively. 
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Note that in the case (3 ^ a 2 , the Hamiltonian Hy cannot be transformed 
into Birkhoff normal form up to order 4 due to resonances. Nevertheless, the 
Hamiltonian truncated at order 4, Hy unc , can be proved to be completely in- 
tegrable. The form of the resonance lattice introduced above suggests that 
Ik + lN-k (1 < k < ^r) are integrals of Hy unc in involution. To find the remain- 
ing commuting integrals we express H a ^(I) in terms of Ik + In-u (1 < k < 4r) 
and a remainder term, 

f-l 

H aJ (I) = ffW(J) + ffW (J) + _L £ rf-47 w _ fe (14) 

fe=i 

where 

f -l 

tf( 2 )(7) - 2 £ sin ^(4 + I N _ k ) + 21 n , 



fc=i 



/3 - a 2 T ^ . fcTr 



N - ^ N 

k=l 

,,2 



1 

7« ^ sin— (4 +7jv-/, ) 



/3 — or . kir , lir 

+ 9/v 2^ sm — sm— (4 + / Ar _ fe )( / ; +^-0, 



and 



2iV ^ TV AT 

i<f=,i<f 



d fe = d fc (a, /3) := -a 2 + (/3 - a 2 ) sin 2 



By ([7]), one has I k lN-k = Jk + Mj* f° r an y l<fc<-y~lso that the remainder 
term Sfc=i ^fc Ik^N-k m (| 14[) can be written as 

2^( E (^(J, 2 + M 2 ) + < ir(j| + M|))+^(j|+M|) 
\i<fe<f 

where the latter term is defined to be if 4^ ^ Z and fc = fc(fc) = £ — k. 
Combined with the expression (JTTJJ) for R a p(J, M) the Hamiltonian Hy unc in 
([15]) then takes the form 

Ht j Unc= NP^_ +fl(J)(J)+fi W (i)+ 1 £ K k (I,J,M) (15) 



Kfe<J 



A? 



where for 1 < fc < 4 
7T fc (7, J, M) = d-(J 2 +Af 2 )+dr(J?+Af|)~2(/3-a 2 ) sin — (J k J- k -M k M- k ) (16) 



■5 



and 

rsrr.rs i -a 2 J 2 N + (/3 - 2a 2 )M 2 N if^ = 0mod4 
K»{J,M) = { x ~ (17) 

4 I otherwise. 

Theorem 1.2. Let N > 4 fee an even integer. Then the truncated FPU Hamil- 
tonian Hy unc given by H3\) is completely integrable. If the expansion of V in 
|[3J) satisfies (a, [3) 7^ (0,0), then the following N — 1 quantities are functionally 
independent integrals in involution: 

(I k + lN-k)i<k<%, (4 + /« + fc)l<fe<«, {K k )i< k <N- 

Remark: In the case (3 = a 2 , the integrals K^{I, J,M) only depend on the 
action variables /, as M 2 + J 2 = IjlN__- for any 1 < j < y. 

In sections [3] and [4] we present a detailed analysis of the geometry of the 
moment map of the integrable system of Theorem 11.21 In particular we show 
that whenever (3 ^ a 2 , then this integrable system exhibits hyperbolic dynamics. 

In section [SJ we will use Theorem 11.11 to show that any FPU chain with 
Dirichlet boundary conditions admits a Birkhoff normal form up to order 4 by 
recognizing such a system as an invariant submanifold of a periodic FPU chain. 
Consider a chain with N' (N' > 3, not necessarily even) moving particles and 
Hamiltonian given by 

at' n' 

n—l n—1 

Assume that its endpoints are fixed, i.e. 

<7o = QN'+i = 0. (18) 

Theorem 1.3. Any FPU chain with N' > 3 moving particles and Dirichlet 
boundary conditions admits a Birkhoff normal form of order 4, i.e. there are 
canonical coordinates (xk,yk)i<k<N' so that Hy takes the form 

{N '\ 1)P2 ' +H° fj (I) + 0(\( X ,y)\ 5 ) 

where I = ,In') is given by £§J), a, (3 are as in $3$, and H„ is of 

the form 



N' , N' 



/ 



(3 — a 



2 



16(7V' + 1)^ V v 32(AT' + 1) 

only if -j-eN 



[3-d 2 
16(jV' + l) 



N' 

^ 4:Sis m JiI m - ^2 s 2 khlN'+i-k ■ (19) 



k=l 

l<l,m<N f 
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Note that in \19jl , the numbers Sk = sin ^ = sin f or an y ^ < k < N' are 

pairwise different. 

Corollary 1.4. Near the equilibrium state, any FPU chain with N' moving 
particles and Dirichlet boundary conditions can be approximated up to order 4 
by an integrable system of N' harmonic oscillators which are coupled at fourth 
order except if (3 = a 2 . 

Denote by p the Hessian of at 7 = 0. Note that ^ is an 

N' x N' matrix which only depends on the parameters a and (3. 

Theorem 1.5. (i) For any given a G M \ {0}, det(Q^ > /3 ) is a polynomial in 
(3 of degree N' and has N' real zeroes (counted with multiplicities). When 
listed in increasing order, the zeroes f3k — Pkipi) satisfy 



Pi < • ■ ■ < (3 r N!+i^ <a 2 < (3 r N>+3^ < ...< (3_ 



N' 



Moreover index(Q® defined as the number of negative eigenvalues of 
Qa,f3> is 9 iven °y 

for p<fa 

index{Q® fj ) = { for /3 r «^+i n < (3 < f3 rN >+^ 

for (3 > f3 N , 




(ii) For a — 0, det(Q^) is a polynomial in (3 of degree N' , and f3 = is the 
only zero o/det(Q^). It has multiplicity N' , and the index of Qq a is 
given by 

for f3 < 



index (Q 0/3 ) = 




for 13 > 



Applications: Theorems 11.31 and 11.51 allow to apply for any given a € K the 

classical KAM theorem (see e.g. 0) near the equilibrium point to the FPU chain 

with Hamiltonian Hy and Dirichlet boundary conditions for a real analytic 

potential V(x) = \x 2 - f x 3 + §x 4 + . . . with (3 e M \ {/?i(a), . . . ,(3 N ,(a)}. 

Moreover, as for any given a € M\{0}, 3 is positive definite for (3 w-i (a) < 

L 2 J 

(3 < [3 n'+i {a), one can apply Nckhoroshev's theorem (see e.g. [TU]) near the 

L 2 J 

equilibrium point to the FPU chain with Hamiltonian Hy for V with such /3's. 
These perturbation results confirm long standing conjectures - see e.g. [3]. 

In subsequent work we plan to compare orbits of the FPU chain with Hamil- 
tonian Hy with the ones of its integrable approximation Hy unc , with the in- 
tention to explain the numerical experiments of Fermi, Pasta, and Ulam. 

Related work: Theorem 11.21 and Theorem 11.31 improve on earlier results of 
Rink in [11] and [13], respectively, where the case a = has been treated. Our 
approach has been shaped by our work on the Toda lattice [6]. The latter one, 
introduced by Toda Q~l] and extensively studied in the sequel, is a special FPU 
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chain which is integrable (i.e. the original Hamiltonian Hy - not only Hy unc - 
is integrable). It turns out that the same canonical transformations which near 
the equilibrium bring the Toda lattice into Birkhoff normal form can be used 
for any FPU chain. If N is odd, this transformation brings Hy into Birkhoff 
normal form up to order 4; we have studied this case in detail in [7]. 

One of the most important open problems in the field of FPU chains is the 
investigation of the dynamics of these chains when the number of particles gets 
larger and larger. It is likely that our results on FPU chains with Dirichlet 
boundary conditions can be used for this purpose. For recent contributions in 
this direction see e.g. [I], [2]. 

Acknowledgement: It is a great pleasure to thank Bob Rink for valuable 
comments on earlier versions of this paper, and Gerda Schacher for her help 
with the graphics. 



2 Proof of Theorem [L2] 

Denote by {•, •} the standard Poisson bracket on M 2Ar ~ 2 . In a straightforward 
way one computes the Poisson brackets between the variables /, M, J,L £ I™" 1 , 
given by © and © (cf. 4 , p. 28): 

Lemma 2.1. The Poisson brackets between the variables I k , Jk, Mk A < k < 

N — 1 ) are given by 

{I h h} = {J h J k } = {M h M k } = 0, (20) 
{J l: Ik} = -Mi(Ski-S k+LN ), (21) 
{M h I k } = Ji(6ki-6 k+ i, N ), (22) 

(23) 

As a consequence, one obtains the following relations between the variables Mk, 
Jk, and Lk, 1 < k < N — 1: 

{M k , Ji} = L[(6 k +i,N - Ski), 
{Jk,U} = M k (6 k +LN - Ski), 
{L k ,Mi\ = Ji(5k+i,N — hi)- 
First note that the list of functions of Theorem 11.21 

(7fc + Jjv-fc)i<fc<^, {h + I%+k)i<k<%> ( K k)i<k<l$.> C 24 ) 

contains N— 1 terms regardless whether ^ is an integer or not. In addition, for 
any 1 < k < the terms I k + lN-k, In /2-k + In /2+k, h+lN/2+k, K k are func- 
tions of the eight variables Xk,y k ,x N / 2 _ k ,y N / 2 _ k ,x N / 2+k ,y N /2+ k ,x N ^k, and 
yN-k, the term In/2 is a function of the two variables xn/2,Vn/2, and, in the 
case N 74 € N, the terms In. + Jm , Kn_ are functions of the four variables 

' 4 4 4 

XN/A,yN/Ai x ZN/i->D3N/A- Hence we obtain a partition of the 2N — 2 variables 
x\, j/i, ... , Xn—i, Dn-1 m t° l t j 1 pairwise disjoint sets of variables, and all 
Poisson brackets between variables of different sets of this partition vanish. 
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Lemma 2.2. The N — 1 functions listed in \2J$ are pairwise in involution. 



Proof. The functions in (|24|) depend on only one of the i_f j + 1 pairwise disjoint 
sets of variables. As the Poisson brackets between terms depending on variables 
of different sets vanish, it remains to check that functions of with the same 
k are in involution with each other. In view of the definitions (HHJ) and (TIT)) 
of Ki and ifjv/4) respectively, and taking into account that (Ik)i<k<N~i are 
pairwise in involution, this amounts to proving that for any 1 < I < 

{JiJf-i ~ MiMk_ u h + Itf-i} = 0, (25) 

{ Ji Jf - M t M § _„ Is _ x + = 0, (26) 

{JiJ,_ t - M,M f _„ I t + I f+l } = 0, (27) 

and 

{ J« - M« , In + I 3N } = 0. (28) 

~4~ ~4" 4 4 v ' 

First we note that by (dU) and flU]) one has for any 1 < I < ^ 

UiJf-i ~ MiMn_ 1 Ji} = -Jn^Mi - M|_, J z (29) 

and 

{JiJ^-i ~ MiMn^I n ^} = Jn_ 1 Mi+ MN_ t Ji. (30) 

Since the right hand sides of (f2"9"| and (|3H|) are invariant under exchanging I and 
Y — I, the same must hold for the left hand sides, and we conclude that 

{JiJn^i - MiMn _;,/«_,} = -J.£_jMj - Mn_ v Ji (31) 

and 

{J/J«_ ; - MjMn_j,7w; + J = Jn_ 1 M [ + M|_,Ji. (32) 

The identities (|25p - ([2T]) now follow from the appropriate combinations of (f2"9"| - 
In the same fashion, one concludes that J28D holds. □ 



Proof of Theorem \1. 6 A In view of Lemma [2. 2 [ it remains to check that the quan- 
tities listed in (|24|) are functionally independent integrals. The independence is 
easy to verify, and the fact that they are conserved quantities follows from the 
formula (fl5|) . showing that Hy unc can be written as a function of them. □ 



3 Foliation by Liouville tori 

In the next two sections we describe the geometry of the moment map of the 
truncated resonant normal form (|13[) for any even periodic FPU chain Hy with 
potential V whose expansion §5§ satisfies (a,/3) ^ (0,0). The case (3 — a 2 is 
special as in this case the normal form (|13[) is the Birkhoff normal form of order 
four of the Toda chain. Its foliation is well known - it is the one of uncoupled 
harmonic oscillators. Hence we will concentrate on the case (3 ^ a 2 only. The 



9 



special case a — has been partially studied by Rink [T2] . Surprisingly, it turns 
out that many of his results continue to hold in the general case. Using the 
notation fc = k(k) = ^ — fc, the integrals of Theorem 11.21 can be grouped as 
follows: 

(CK k ,CK k , L k , K k ) 1 < k< r^, 7w, ^Njifw, (33) 
where for 1 < k < ^ 

'■= Ik + iN-k, Lk '■= Ik — !%■ 

(Here we used that L k = (I k + Iq.+ k ) ~ (1% -fc + ^-f +fc) ^ s tnc difference of two 
integrals listed in Theorem 1 1.21 ) 

Using the assumption (3 — a 2 ^ 0, we rewrite the integrals (i^/) 1 </<» as 
follows. Introduce the bifurcation parameter 

1 = 1 {a,[3):=-^- (34) 
u z — (3 

and note that 

d- = -a 2 + (/3 - a 2 )4 = {13- a 2 )( 7 + 4)- 
For any 1 < fc < -j, one has 

K k = -d-(J 2 + M 2 ) - d- k {J 2 k +Ml) + 2{f3-a 2 )s 2k (J k J k -M k M k ), 
whereas for k = -j, 

K k = a 2 , J 2 -{13- 2a 2 )M 2 = {a 2 - (3) {jj 2 + (1 + ^M 2 ) . 

In the sequel, we will for simplicity omit the factor S2 k (a 2 — 13) in K k , since it 
does not influence the geometry of the level sets of the integrals {K k ) l<k< N_ . 

Each of the l|j+1 groups of integrals listed in (j3"3"|) depends only on a subset 
of the variables {(x k , y k )i< k <N-i}- These subsets form a disjoint partition of 
{(%k, Uk)i<k<N-i}- More precisely, the following result holds. 

Proposition 3.1. The phase space T*M. N ~ 1 of the truncated resonant normal 
form Hy unc given by \1S\) is the direct sum of invariant symplectic subspaces 

T*R N ~ 1 = ?fc, 

0<fc<f 

where 

y k = {(x j ,y j )i< j < N -i eri"- 1 ^ = yj =ovi^ {fc,iv-fc,fc,iv-fc}}. 

The foliation of T*]^ -1 by level sets of the integrals t33\) is the Cartesian 
product of the foliations of the y k . 
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We now analyze the foliations of To, fn, and for < fc < -j- separately. 
Note that and CPfc (0 < k < -j) can be canonically identified with T*R and 
T*M 4 , respectively, and Tw with T*R 2 (if f e N) or {0} (if f $ N). 

Foliation of ?q One easily sees that Jw. foliates T*]R by circles, centered at the 
origin. 

Foliation of Tw /or ^ e Z Let us study the geometry of the moment map 
M : T*R 2 — > R 2 defined by the integrable system with commuting integrals H = 
H n and K = K n . It is convenient to introduce the following notation. Denote 

4 4 

the standard coordinates of T*R 2 by (x,y) = (xi, x 2 , J/i, 2/2) and introduce the 
action variables Ij = \{x 2 (j = 1,2), as well as the Hopf variables M, J, L 

given as in ©, 

(M,J,L) = ^{x!y 2 - x 2 yi,x 1 x 2 + yiy 2 ,Ii - h)- 

Then the moment map M = (H, K) takes the form 

H=-(h+I 2 ) and X= (1 + 7 )M 2 + 7 J 2 . 

As already remarked in ([7|) one has 

M 2 + J 2 + L 2 = H 2 . 

Further, we may replace K by K-y given by 

f (1 + 7 )M 2 + 7 J 2 7 £ {0, 1} 
if 7 := <^ M 7 = 

I J 7=1 

First observe that the origin (x, y) — (0, 0) of T*M. 2 is the only critical point of 
M with rank^^^M) = 0. Moreover, 

M- 1 {(0,0)} = {(0,0)}. 

The critical points (x, y) £ T*R 2 \ {(0, 0)} with rank(d( Xj j,)M) = 1 are analyzed 
by symplectic reduction via the Hamiltonian vector field of H . On the sphere 
§^ = {H = p 2 /A} of radius p > in T*M 2 define the Hopf map 

3:S 3 p ^S 2 r ,(x,y)^(M,J,L) 

where r = V M 2 + J 2 + L 2 \§3 = H\g3 = B-,. The fibers of J are circles obtained 
by the S 1 -action of H. The reduced system is then given by (§ 2 , X 7 ) where X 1 
denotes the reduced Hamiltonian vector field induced by K 1 . To compute X 7 , 
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note that the equations of motion in the reduced system corresponding to the 
Hamiltonian K 7 are given by 



(35) 





Indeed, following the procedure of reduction in section 1.5 of [3], formula ([35 
follows from 



i} (1 < 3 < 3) 



i=i 



and the commutation relations of the variables (twi, 102,103) = (M, J, L) given 
by Lemma [2TTI We then obtain 

(-27JZ, 2(1 + i)ML, -2MJ) 7 i {-1, 0} 

(0,L,-J) 7-0 
(-L,0,M) 7 = -l 



It turns out that the foliation of E 2 . by level sets of K 1 depends on the bifurcation 
parameter 7. If 7 = 0, then (±r, 0, 0) are the only two fixed points of X . They 
are both elliptic and the level sets of K in E> 2 \ {(±r, 0, 0)} are circles. Similarly, 
for 7 = — 1, (0, ±r, 0) are the only two fixed points of X_\. They are both elliptic 
and the level sets of K-\ in E> 2 \ {(0, ±r, 0)} are circles. Now let us consider the 
case 7 € M \ {0, —1}. Then X 1 admits six fixed points, 

(±r,0,0), (0,±r,0), (0,0, ±r) 

where two of them are hyperbolic and the remaining four elliptic. Note that 
the corresponding values of K 1 are (1 + j)r 2 , jr 2 , and 0, respectively, and that 
the two hyperbolic fixed points are contained in the same connected component 
of the inverse image of K 1 in § 2 . This component consists of two great circles 
where each of the four half circles is a heteroclinic Xy-orbit connecting the two 
hyperbolic fixed points. 





hyperbolic fixed points 


critical value 


7 < -1 


(±r,0,0) 


{l+j)r 2 


-1< 7 < 


(0,0, ±r) 





7 > 


(0,±r,0) 


>\ 



(36) 



Let us verify the claimed classification of the two fixed points (0,0, er) with 
e E {±}- The other four fixed points are treated in a similar fashion. Near 
(0, 0, er) we choose M, J as coordinates of S> 2 . The equations of motion induced 
by on E> 2 in these coordinates read 



M 
j 



-2-fJ^r 2 - M 2 - J 2 , 



e2(l + i)M^r 2 - M 2 - J 2 . 



12 4 INTEGRABLE SYSTEM CK K ,% R ,G K ,K K ) 

If linearized at (0, 0, er) the corresponding linear system is given by the 2x2- 
matrix eA where 

A = 2r(,° 7' 
+ j 

The eigenvalues of A are given by 



A 2 +4 7 (1 + 7 )r 2 = or A x , 2 = ±2ry/-j(l + t)< 

As -7(1 + 7) > iff -1 < 7 < it follows that Ai,2 are in R \ {0} and hence 
that (0, 0, ±r) are both hyperbolic fixed points for —1 < 7 < whereas they are 
both elliptic if 7 < —1 or 7 > 0. For — 1 < 7 < 0, the inverse image i4T~ 1 ({0}) 
is given by 

KyH{0}) = i( M > J > L )\( 1 + 7)M 2 + 7 J 2 = 0; M 2 + J 2 + L 2 = r 2 } 

1/2 



{(M,J,L)\M = ± 
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J; M + J + L = r 2 }, 



1 + 7 

whereas for 7 < -1 or 7 > 0, J fif- 1 ({0}) = {(0,0, ±r)}. 

Foliation of Tf. for < fc < x The foliation of for < k < is more 
complicated than in the case k = ^ . We describe it in detail in the subsequent 
section. 



4 Integrable system (Jik, ^Ki, Gk, Kj~) 

In this section we present a detailed study of the geometry of the moment map 
M : T*R 4 — * R 4 defined by the integrable system 3^ with commuting integrals 
JCfc, JCj,, Gk, Kk for 1 < k < As in section [3] we restrict to FPU chains with 
potential V whose expansion (J3j> satisfies (3 ^ a 2 . We show that in this case 
the vector field induced by Kk exhibits hyperbolic dynamics. It is convenient to 
introduce the following notation. Denote the standard coordinates of T*R 4 by 
(x,y) with x = (iEi)i<i<4 and y — (j/i)i<i<4, and introduce the action variables 
Ij = \{x 2 + y 2 ) (1 < j < 4), as well as the Hopf variables (Mi, Ji, ij)i<i<2 
given by 

(Mi,Ji,Li) = -(xij/2- #22/1, £12:2 + 2/12/2,-^1 —h), 

(M 2 , J 2 ,L 2 ) = -(^32/4 - ^42/3, ^3^4 + 2/32/4, ^3 - h)- 

By Lemma [2~T1 the Poisson brackets between the variables (Mi, Ji, Li)i<i<2 arc 
given by 

{Mi,Ji} = -Li, {Ji,Li} = -M h {L h Mi} = -Ji 



whereas all other brackets vanish. 
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The moment map M then takes the form 

M:T*R 4 ^R 4 , {x,y)^(H 1 ,H 2 ,G,K^) 

where 

H 1 = ^(I 1 +I 2 )\ H 2 = i(/ 3 +7 4 ); G = L 1 - 
and where K y is a scalar multiple of Kk, given by 
2 i 

i=l 

with the coefficients di, 7 , d2, 7 defined by 



di 7 — 



7 + 4 , 

S2fc 



C?2,7 — 



7 + 4 

S2/c 



and Sfc = sin^£, = cos^-. (The definition of the integral G above differs 
from the one given earlier by a multiple of the integral Hi — H 2 as I\ — 13 = 
L\ — L2 + Hi — i?2-) 

First note that the origin (0, 0) in T*R 4 is the only critical point of M with 
rank(cDVt) = 0. Moreover, M -1 {( >°)} = {(0,0)}. Next observe that when 
restricted to T*R 2 x {0}, one has G = \{h - h) and K 7 = di, y (Hf - L 2 ), 
hence they are functions of Ii, I2 alone and M|t»r2 x { } may be replaced by the 
map (x,y) 1— ► (ii,/2,0, 0). The geometry of the latter map is the one of two 
uncoupled harmonic oscillators. The subspace {0} x T*R 2 is treated similarly. 
It remains to study the restriction of JVC to T*R 4 \((T*R 2 x {0})U({0} x T*R 2 )). 
The Hamiltonian vector fields of Hi and H2 induce a torus action on T*R 2 . The 
corresponding symplectic reduction is given by the product of two Hopf maps, 

?:§3 x§ 3 _^ § 2 x§ 2 (^^^(Mi.Ji.LOi 



)l<i<2 



where for i = 1,2, 



= {Hi = pf/4:} is a sphere in T*R 2 and r t = pf/4. = 
Hi\§3 . The fibers of 5F are 2-dimensional tori, ob- 
The reduced system is then given 



y/Mf + J-f + L^ 

tained by the S 1 x § 1 -action of Hi x H 2 

by (S~ x S? , V, Xy), where Y and Xy denote the reduced Hamiltonian vector 
fields induced by G and K 7 , respectively. To compute Y and Xy, note that the 
equations of motion in the reduced system, corresponding to a Hamiltonian H , 
are given by 






1,2 



(37) 



see section [3] for details. We then obtain 

/ Ji \ 
-Mi 



Y = 




-h 
M 2 

V j 



X<y 



/ ( J% - di, 7 Ji)Li \ 
{di,jMi + M 2 )Li 
-(M 1 J 2 + M 2 Ji) 
(Ji — d<i a J<i)Li 
(d2, 7 M 2 + Mi)L 2 
\-(M 1 J 2 + M 2 Ji) J 



(38) 
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Further introduce the reduced moment map 

M 7 :S 2 n x S£, -^R 2 ,(M 4 , J l ,L l ) 1 < % <2 >-» (G,K y ). 

We now study the critical points of M 7 with rank (dM 7 ) = 0, i.e. points of 
§ri x ^r 2 wm ch are fixed points of both, Y and X 1 . From the expressions for 
Y and X 7 derived above, one easily sees that there are only four such critical 
points, 

{Mi, J h L 4 )i< 4 < 2 = e(0, 0, n, 0, 0, ±r 2 ), 

where e € {±}- The value of the critical point e(0, 0, ri, 0, 0, — r 2 ) by M 7 is 
(e(ri + r 2 ),0) and 

M^Hfcfa +r 2 ),0)} = {£(0,0,n,0,0,-r 2 )}. 

Computing the Jacobian of the reduced vector field X 1 at the critical points 
one sees that they are elliptic fixed points of X 1 . The values of the other two 
critical points e(0, 0, r\, 0, 0, r 2 ) by M 7 are {e(ri — r 2 ),0). The inverse image 
of (e(ri — r 2 ),0) might have several connected components, depending on the 
values of 7 and the additional bifurcation parameter 

r:=^>0. 

Our main results concerning the critical points e(0, 0, r\, 0, 0, r 2 ) are collected 
in the folloing theorem. 

Theorem 4.1. Assume that 1 < fc < f , < r < 1, ee {±}, and 7 g K. The 
critical point e(0, 0, r±, 0, 0, r 2 ) 0/ M 7 is a hyperbolic /ized point of the vector 
field X 1 if and only if 



( 7 +4)vf+(7+4)_L 



< 2s 2k . (39) 

Otherwise it is an elliptic fixed point of X^. If 139\) is satisfied, the stable and un- 
stable manifolds o/e(0, 0, 7*1, 0, 0, r 2 ) ooi/i /ia?;e dimension two. In the case r < 1, 
the connected component o/M~ 1 {e(ri — r 2 ,0)} containing e(0, 0, ri, 0, 0, r 2 ) is 
a 2-dimensional torus pinched at e(0, 0, ri, 0, 0, r 2 ) and consists of homoclinic 
X 1 -orbits. In the case r = 1, M^KOjO)} is a 2-dimensional torus pinched at 
the two points ±(0, 0, 7*1, 0, 0, ri), and M~ 1 {(0, 0)}\ {±(0, 0, r 1; 0, 0, r^)} consists 
of heteroclinic X^-orbits. 

To prove Theorem 14.11 we separately treat for any given 1 < k < three 
subsets of the domain of the parameters 7 and r. The results for these three 
cases are stated in detail in Propositions 14.21 14.31 and 14.41 below. 

Note that the inverse image M~ {(e(ri — r 2 ), 0)} is invariant under the action 
of the vector field Y. The orbits of this action can be easily described, 

{(R-^Mt, Ji),L 1 ,^(M 2 ,J 2 ) ) L 2 )||0| <tt}, (40) 
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(c) k = 10 (d) k = 24 



Figure 1: Subsets of parameters (7, r) with hyperbolic dynamics of X~, (iV= 100) 

where R<f,(u, v) denotes the image of (u, v) e K 2 of the rotation by the angle 
<j> in counterclockwise orientation. Hence given L\ and Z/ 2 with |Lj| < rj for 
i = 1,2 there exists a unique point (Mj, Ji, £j)i<j<2 on such an orbit with the 
property that 

(Mi, Ji) = (jrf - L\, O^j and (M 2 , J 2 ) = ^/r| - i| (cosa,sina) 
for some < a < 27r. We denote the corresponding F-orbit by H(L\, L 2 , a), i.e. 

£(ii, i 2 , a) = { (R-^rf-L 2 v 0), Li, R a+<f> {^ rl~ L% 0), L 2 ) 1 10| < tt}. (41) 

As G and i"T 7 commute, K y is invariant along any orbit of the vector field Y 
and we conclude that 



if 7 ((M, Ji, L I )i< l < 2 ) = - Y,d in (rf - L\) + yJ(rl-Ll)(r*-Ll)casa. (42) 

Let us now determine M~ 1 {(e(ri — r 2 ), 0)} for e = 1 and r < 1. (The case where 
r > 1 and/or £ = — 1 is treated in a similar fashion.) Let (Mj, Jj, £,)i<,<2 be 
an element of M" H^i - r 2 ,0)} \ {(0, 0, n, 0, 0, r 2 )}. Then r 2 - L 2 = n - Lj 
and if 7 ((Mj, L i )i< i < 2 ) = 0. First note that L\ < r\. Indeed, if L\ = i\, 
then L 2 = r 2 and (Mj, Jj) = (0, 0) for i = 1,2, contradicting our assumption on 
the point considered. Hence in the expression for K 7 displayed above we can 
factor out 7*1 — L\ and the equation K 7 — reads 

= di, 7 (ri + Li) + d 2n (r 2 + L 2 ) + 2y/(n + Li)(r 2 + L 2 ) cosa. (43) 
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Next let us consider the case where L\ — —r\. Then (Mi, J\,L\) = — (0, 0, ri) 
and (|43|) reads (7 + c^)(r2 + £2) = 0. Hence either L 2 = —r 2 or 7 = — c^. In 
the case L 2 = — r 2 it follows from r 2 — £2 = r i — -^1 = 2r*i that r 2 = ri and 
L 2 = L\. As a consequence (M 2 , J 2 , L 2 ) = — (0, 0, ri). On the other hand, if 
7 = — c\ and r\ < r 2 , then 

-r 2 < r 2 - 2ri = L 2 < r 2 and M 2 2 + Jf = r\ - L\ = 4ri(r 2 - r x ). 

If ii 7^ — ri (and hence L2 + 1 as ri < r%) we set for — r\ < L\ < r\ 

Q(Li) := { V a^+itin if r i < r 2 (44) 
[ 1 if ri = r 2 

Then < Q < \/r and for r < 1, Q is monotonically increasing on — n < Li < 
n . After division by ^- ^/ (ri + L\){r 2 + L 2 ) the equation |43|) reads 

(7 + 4)Q + (7 + c'-)^ + 2s 2fe cos a = 0. (45) 

To investigate the solutions of (|4"5|) we distinguish between three cases: r = 1, 
[7 + c| = and r < 1] , and [7 + c 2 . 7^ and r < 1] . 

Let us first treat the case r = 1. Then Q = 1 and equation (j4"5|) takes the 
form 

27+ 1 = -2s 2fc cosa, (46) 

which is independent of L\. 

Proposition 4.2. Let 7 € R 6e arbitrary and assume 1 < fc < ^ and r = 1. 

TTien i/ie following statements hold: 

(1) // 1 2 7 +l| > 2s 2fc , iften 

M^HCO, 0)} = {e(0, 0, n, 0, 0, n)\e = +}, 

and ±(0, 0, ri, 0, 0, n) are ooi/i elliptic /ixed points of the vector field X 1 . 

(ii) If 1 27 + 1| < 2s 2 fc, i/ien M~ 1 {(0, 0)} = JST a U N-a, where a is the unique 
angle satisfying < a < ir and 2j + 1 = —2s 2 k cos a, and where for any 

-IT < (3 < IT 

!N> = |J L{L u L 2: f3) 

ill<'\L 
L 2 = L 1 

with L(Li, L 2 , (3) given by \4ty . Both points, ±(0, 0, ri, 0, 0, ri), are hy- 
perbolic /raed points of X l7 and their stable and unstable manifolds have 
each dimension two. The set N Q \ {±(0, 0, r%, 0, 0, r±)} consists of het- 
eroclinic X^-orbits from (0, 0, n, 0, 0, r%) to — (0, 0, r\, 0, 0, r\), whereas 
N_ Q \ {±(0, 0, r\, 0, 0, ri)} consists of heteroclinic X~ i -orbits with opposite 
direction. Topologically, M~ 1 {(0, 0)} is a 2-dimensional torus, pinched at 
each of the two fixed points ±(0,0, 7*1, 0,0, n). 
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(m) 7/27 + 1 = -2s 2k , then a = and M~ 1 {(0, 0)} = X , w/iereas z/2 7 +l = 
2s2k, then a = ir and M~ 1 {(0, 0)} = N^. In both cases, ±(0, 0, rx, 0, 0, rx) 
are elliptic /Kced points of X 1 . On N U N^, any -orbit is periodic and 
coincides with the corresponding Y -orbit at least up to orientation. 

Proof, (i) By a straightforward computation one shows that under the given 
assumptions, both points ±(0, 0, r%, 0, 0, r±) are elliptic fixed points of X 1 . In 
view of equation (|46|) . item (i) then easily follows. 

(ii) By the discussion preceding Proposition 14.21 it follows that the inverse 
image M~ 1 {(0,0)} is given as claimed. Again by a straightforward compu- 
tation one shows that both fixed points ±(0, 0, rx, 0, 0, n) are hyperbolic. To 
see that N Q \ {±(0,0, ri, 0,0, ri)} consists of heteroclinic orbits of the vector 
field Xy, consider the third component (Xy)3 of X 1 (cf (|38|). Any element 
(Mj, Jj, Z/j)i<i<2 in N a is of the form 

(M x , Ji) = yJrl-Ll(cos <f>, - sin 0), (M 2 , J 2 ) = \J rl-L\{cos(a+4>), sin(a+0)). 

(47) 

Thus 

(X 7 ) 3 = -{M 1 J 2 + M 2 J 1 ) 

= -(r\ - L\) (cos0 sin(a + 0) - cos(a + 4>) sin</>). (48) 

Hence (Xy) 3 = — (rf — L?)sina < for any point in N a \ {±(0, 0, rx, 0, 0, n)}. 
As the last component of X 7 coincides with the third one, it follows that any 
Xy-orbit on N Q \ {±(0, 0, ri, 0, 0, r%)} originates from (0, 0, r%, 0, 0, rx) and ends 
in —(0, 0, rx, 0, 0, ri). The orbits on N_ Q are analyzed in a similar way. 

(iii) Clearly, if 2j + 1 = — 2s2k, one has M~ 1 {(0, 0)} = No and one verifies 
in a straightforward way that ±(0, 0, rx, 0, 0, r\) are elliptic fixed points of X 7 . 
According to l|48p. the third component {X^)^ of X 1 vanishes identically on No. 
Further, 2 7 + 1 = —2s 2 k implies that 1 + d 2 . 7 = —1 — di. 7 . In view of l|47p it 
then follows that 

X 7 = (1 + d 2n )L 2 ■ Y. 

The claimed statements for the case 2 7 + 1 = 2s 2 k are proved in a similar 
fashion. □ 

Next we consider the case where r < 1 and 7 + c\ = 0. Then 7 + s\ = —c 2 k 
and hence 

di 7 = and d 2 7 = 0. 

S2fe 



Thus equation (|45|) takes the form 

c 2 kQ(Lx) = 2s 2k cosa. (49) 
Note that < c 2fc < 1 as 1 < k < f . 

Proposition 4.3. Assume that 1 < k < < r < 1, and 7 + c\ = 0. TTiera 
the following statements hold: 
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(i) If yfr > 2s2k/c2k, then the connected component of M.~ 1 {(ri — r2,0)} 
containing the critical point (0, 0, r\, 0, 0, r 2 ) consists of this point alone. 
It is an elliptic fixed point of X 7 . 

(ii) If y/r < 2s 2 fc/c 2 fc, then 

M- 1 {{r 1 -r 2 ,0)}= \J C J (L 1 ,L 2 ,a Ll )UL(L 1 ,L 2 ,-a Ll ) 



[£ii<n 

^2 = I 'l + r 2- T "l 



where for any \Li\ < n, (Xl 1 is the unique angle satisfying 

7T 

C2kQ{Li) — 2s 2 fc cosa Ll and < a Ll < — 

and Q(Li) denotes the function defined by continuously extended to 

the closed interval [—L\,L-\\. Furthermore, the connected component of 
M~ 1 {(ri — ?"2,0)} containing (0, 0, r±, 0, 0, r 2 ) consists of homoclinic X~ ( - 
orbits which originate and end in (0, 0, r±, 0, 0, r 2 ). Topologically, it is a 
2-dimensional torus, pinched at (0, 0, n, 0, 0, r 2 ). 

If y/r < 2s 2 fc/c 2 fc 7 £/ien (0, 0, ri, 0, 0, r 2 ) is a hyperbolic fixed point of X^ 
and its stable and unstable manifold have each dimension two. If y/r — 
2s 2 fc/c 2 fe, (0, 0, ri, 0, 0, r 2 ) is an elliptic fixed point of X 7 . 

Proof, (i) By a straightforward computation one shows that under the given 
assumptions, (0, 0, n, 0, 0, r 2 ) is an elliptic fixed point of A 7 . In view of equation 
and the discussion of the case L\ — ±n item (i) then follows easily, 
(ii) Again by the discussion preceding Proposition 14.21 it follows that the 
inverse image M~ 1 {(ri — r 2 , 0)} is given as claimed. Again by a straightforward 
computation one sees that (0, 0, n, 0, 0, r 2 ) is a hyperbolic fixed point of X 1 if 
y/r < 2s2k/c2k, and an elliptic one if y/r = 2s 2 k/c 2 k- Next consider a point 
{MuJ^L-^iKi in M- 1 {(ri-r 2 ,0)} with 



(Mi, JO = \jr\-L\ (cos^, 



— sin ( 



(M 2 , J 2 ) = \jr\-L\ (cos(a Ll + 4>),sm(a Ll + 0)) 
where \L\ \ < r\. Then the third component of X 1 is given by (cf (08])) 



(A 7 ) 3 = -yjrf -L\\jr\-L\ sina Ll . 

Hence (A 7 ) 3 = (A 7 ) 6 < 0. It follows that the A 7 -orbit passing through such 
a point originates at (0, 0, r\, 0, 0, r 2 ) and then reaches a point of the form 
(0,0,-ri,M 2 , J 2 ,L 2 ) with 

L 2 = r 2 -2ri > -r 2 and (M 2 ,J 2 ) = \jr\- L|(cos(7T + 0), sin(7T+0)). (50) 
At this point the vector field A 7 is given by 



(-riJa.-nMa.O.O.O^). 
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Note that this vector does not vanish as M| + Jf = r\ — L\ > 0. Similarly, at 
a point (Mi, Ji, Lj)i<j< 2 in M^ 1 {(r i - r 2 , 0)} satisfying 



(Mi, Ji) = ^r\-L\ (cos0, -sin</>), 
(M 2 , J 2 ) = \jr\-L\ (cos(-a Ll + 0), sin(-a Ll + </>)) 
and |Li| < ri one has 



(Xy) 3 = \Jr\-L\^r\-L\ sin a Ll ■ 

Hence (^7)3 = (^7)6 > 0. It follows that the X 7 -orbit passing through such 
a point ends up at (0, 0, 77, 0, 0, r 2 ) and passes through a point of the form 
(0, 0, — n, M 2 , J 2 , -^2) w hh (M 2 , J 2 ,L 2 ) as in ([50)) . We then conclude that the 
connected component of M~ 1 {(ri — r 2 ,0)} containing (0, 0, 77, 0, 0, r 2 ) consists 
of homoclinic X 7 -orbits originating and ending at (0, 0, 77, 0, 0, r 2 ). □ 

Finally let us treat the case r < 1 and 7 + 7^ 0. Denote by Q(Li) the 
function defined by (|44[) . extended continuously to the closed interval [— 7*1, 7*1]. 
Further introduce the function 

/ : (0, VF) - R, q » (7 + 4)g + (7 + C 2)I. (51) 

q 

Note that lim 9 \^ 1/(9)1 =00. 

Proposition 4.4. Assume that 1 < fc < < r < 1, and 7 + c| 7^ 0. T/ien 
£/ie following statements hold: 

0) If \f{y/r)\ > 2s 2 fc, i/ien t/ie connected component of Jd~ 1 {(ri — r 2 ,0)} 
containing the critical point (0, 0, n, 0, 0, r 2 ) consists of this point alone. 
It is an elliptic fixed point of X 7 . 

(ii) If\f(s/r)\ < 2s 2 fe, then there exists —77 < l~ t , r < 7*1 so that the connected 
component o/M~ 1 {(ri — r 2 ,0)} containing (0, 0, 77, 0, 0, r 2 ) is given by 

[J L(L 1 ,L 2 ,a Ll ) U Z(L 1 ,L 2 ,~a Ll ) 

i 2 = - L l+ T '2- 1 ~l 

where for any l 7 _ r < L\ < 77, «ij is £/ie unique angle satisfying < ai^ < 
7r and 

f(Q(Li)) = -2s 2fe cos(a Ll ). 

TVie point (0, 0, 77, 0, 0, r 2 ) is a hyperbolic /ixed point of X 7 and its stable 
and unstable manifold each have dimension two. The connected component 
o/M~ {(77 — r 2 , 0)} containing (0, 0, 77, 0, 0, r 2 ) consists of homoclinic X-y- 
orbits which originate and end in (0, 0, 77, 0, 0, r 2 ). Topologically, it is a 
2-dimensional torus, pinched at (0, 0, 7*1, 0, 0, r 2 ). 
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Proof, (i) By a straightforward computation one shows that under the given 
assumptions, (0, 0, r\, 0, 0, r 2 ) is a (possibly degenerate) elliptic fixed point of 
Xy. We have already seen that under the given assumption M~ 1 {(r 1 — r 2 , 0)} \ 
{(0, 0, n, 0, 0, r 2 )} consists of the set of points (Mi, Jj, L^)x<%<2 satisfying L 2 = 
^i + r 2 ~ r i an d PS")) - Note that equation (j4"5")) admits a solution a for Q = 
y/r iff |/(V^)I < s 2k- In the case |/(-\A\)| > 2s 2 fc it follows immediately that 
M- 1 {(ri-r 2 ,0)} = {(0,0,ri,0,0,r 2 )}. If |/(V?)I = 2s 2fc , then an analysis of 
the graph of / near (y/r, f(y/r)) leads to the claimed result. 

(ii) As lim g \^ |/(<z)| = oo it follows that there exists —r\ < Z 7 , r < r\ so that 
the interval [^ n ri] is a connected component of (/ o <5) _1 ([— 2s 2 k, 2s 2 fc])- It 
follows that for any £ 7ir < L\ < r\ there exists a unique angle < < n so 
that 

}{Q{Li)) = -2s 2k cos(a Ll ). 

The connected component of the preimage M~ 1 {(ri — r 2 ,0)} containing the 
point (0, 0, ri, 0, 0, r 2 ) is then given as claimed. Again by a straightforward 
computation one sees that (0, 0, ri, 0, 0, r 2 ) is a hyperbolic fixed point of X 7 . 
One then can argue as in the proof of item (ii) of Proposition 14.31 to show the 
remaining claims. □ 

Proof of Theorem \4-l\ Theorem 14.11 follows from Propositions 14.21 - 14.41 □ 

It remains to study the critical points of M 7 with rank dM 7 = 1, i.e. points 
of (S 2 x x § 2 2 ) \ {±(0,0,ri,0,0,±r 2 )} where the vector fields Y and A 7 are 
collinear. In view of the formulas (f38|) for Y and Xy, points (Mi, Ji, Li) € 5, 2 . , 
i = 1,2, of this type have the property that the determinant of any 2x2- 
submatrix of the 2 x 4-matrix formed by Y and X 7 vanishes. It leads to the 



following system of equations: 

M 1 J 2 + M 2 J 1 = 0, (52) 

J\L 2 + L x J| - Ji J 2 {dx tl Lx + d 2n L 2 ) = 0, (53) 

M\L 2 + UM\ + M\M 2 (di a Li + d 2>1 L 2 ) = 0, (54) 

Ml JiL 2 - L X M 2 J 2 + JiM 2 (di, 7 Li + d 2n L 2 ) = 0. (55) 



Theorem 4.5. Assume that l<k<^r,0<r<l, and 7 e 1. If a point 
(Mi, Ji, Lj)i<i< 2 € §^1 x ^r 2 \ {^(O: 0, ri, 0, 0, ±r 2 )} is a critical point of M 7 
wii/i ronA; dM 7 = 1 then (M 2 , L 2 ) — X(Mx, — Ji) /or some Aet, and 

(r 1 -i?) 2 L|+(rl-^) 2 i?+2(rf-L?)(r|-Ll)(2i 1 i 2 -(d 1 , 7 i 1 +d 2 , 7 L 2 ) 2 ) = 0. 

Given any point (Mi, Ji,L\) € § 2 \{±(0,0, ri)} i/iere exist at most eight points 
(M 2 , J 2 , L 2 ) 6 S 2 2 \ {±(0,0, r 2 )} smc/i t/iat (Mi, Ji, Li)i<i< 2 is a critical point 
of M 7 wit/i ranfc dM 7 = 1 . 
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(a) k = 1 (b) k = 3 (c) fc = 5 

Figure 2: Sets of solutions (Z x , Z 2 ) of ([M]) for N = 24, r = 1, 7 = -0.5 

Proof. First assume that ii G {±^i}- Then Ji = Mi = 0. Hence (|52")l is 
automatically satisfied and equations ([53]) and ([54]) read J 2 = and M 2 = 0, 
respectively. As a consequence, (Mi, J\,L\) = (0,0, ±ri) and (M 2 , J 2 ,L 2 ) = 
(0,0, ±r 2 ). In view of Theorem 14.11 we thus may assume that \Li\ < r±. Then 
(Mi, Ji) 7^ (0,0). Hence the first equation (152"]) says that there exists A e R 
such that 

(M 2 ,J 2 ) = A(Mi,-Ji). (56) 
The conditions (Mi, Ji, Li) € S;!., i = 1, 2 then imply that A satisfies 

A 2 = 3-4- (57) 
r i _ ^1 

Substituting $MB into (j5U]) - (p)5]) one sees, again using (Mi, Ji) ^ (0,0), that 
(|5T?]) - (l5"oT) is equivalent to 

L-2 + X 2 Li + A(di, 7 Li +d 2 , 7 L 2 ) = 0, (58) 

or, taking squares, {L 2 + A 2 Li) 2 — A 2 (di :7 Li + d 2:7 L 2 ) 2 = 0. Using ([ST]) , the 
latter equation reads 

(r 2 -L 2 ) 2 Ll + (r 2 2 -Ll) 2 L 2 +^^^ 

or, after dividing by r\r\ one gets, using the bifurcation parameter r and the 
normed variables U := Lijvi £ (0, 1) (i — 1, 2), 



r 2 (l-/ 2 ) 2 ; 2 + (l-/ 2 ) 2 ? 2 + 2r(l-/ 2 )(l-; 2 )(2;i/ 2 -(v^di, 7 Zi+^d 2 ,^2) 2 ) = 0. (59) 



Note that for given r and < £1 < 1, the left hand side of ([ST)]) is a polynomial 
in l 2 of degree four, (l\ + d 2 7 (l - + 0(Zg). 

Summarizing, we have shown that for any given point (Mi, Ji,Zi) G § 2 t \ 
{(0, 0, ±ri)}, there exist at most 8 points (M 2 , J 2l L 2 ) € § 2 2 \ {(0, 0, ±r 2 )} such 
that Y and A 7 are collinear. Indeed for any (Mi,Ji,Li) S S 2 ^ \ {(0,0, ±ri)} 
a solution (M 2 , J 2 , L 2 ) e § 2 2 \ {(0, 0, ±r 2 )} of ^-fijft is given by (M 2 , J 2 ) = 
A (Mi, Ji) and L 2 — r 2 Z 2 with A and I2 satisfying ([57]) . respectively ([59]) . As a 
consequence the set of solutions of (l52l) - (j55]) is an algebraic subset of § 2 X x § 2 2 
of dimension at most two. □ 
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4 INTEGRABLE SYSTEM ( "K K , % R , G K , K K ) 




In order to analyze the critical points of M 7 with rank(<iM 7 ) = 1, we perform 
another symplectic reduction. First we pass to the orbit space of the flow of Y 
on §J? x S^ 2 and then to the level sets of G. 

In view of ([40]) , the Y-flow is an § 1 -action. Note that besides L\ and L 2 , 
the quantities a and r are invariant under this S^action, 

a := M X M 2 - J\ J 2 , r := Mi J 2 + M 2 Ji. 

They are related by the identity 

- 2 +T 2 =nK 2 -^ 2 ). (60) 

»=i 

Define 

J (3) : x ^ k 4 

(Ml, J 1 ,L 1 ,M 2 , J 2 ,L 2 ) ^ (L X ,L 2 , a, r) 

where S£ ;= §2^ y {(0,0, ±r 4 )}. Let r denote the image of J( 3 l For any 
element (Li, L 2 , s, t) <= O r we have 

2 

s 2 +t 2 = J](rf - Lj) and \L t \ < r< (1 = 1, 2). (61) 

The fibers of are the orbits of the Y-action on x , i.e. O r coincides 
with the orbit space of the Y-action on S 2 x S^ 2 . As a consequence, any function 
on S 2 , x S 2 2 which Poisson commutes with G factors through r . 

In particular, K y and G factor through r . In fact, K~ i and G, when ex- 
pressed in the variables L%, L 2 , a, r, are polynomials, given by 

2 1 

z=l 

G = Li - L a . (63) 

By reducing the system (G,K 7 ) by the Y-action one obtains a family of inte- 
grable systems with one degree of freedom parametrized by the value c of G. 
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Denote by X 1 ^ c the Hamiltonian vector field induced by K 1 . The fixed points 
of Xy_ c can then be characterized in terms of the bifurcation parameters 7, r, 
and k. 

Note that by ([52]) . the rank-l-points of the reduced moment map M 7 satisfy 
r = 0, and by [[56 ]) - ([57 ]) . a 2 = {r\ - L 2 )(r 2 ~ L\). Hence the image of the set of 
the rank-l-points by 3^ is an algebraic subset of O r of dimension at most one 
- see (l59l). 



By (|60p . a and r are located on a circle of radius y/ (r 2 — L\)(r\ — L|), 



(a, r) = y (r 2 - L\){r\ - L|) (cos 0, sin0), (64) 

where G R/27rZ. The phase spaces, reduced by the Y-action, are now obtained 
by taking subsets of O r corresponding to level sets of G, i.e. by replacing L 2 
by L\ — c, where c is the value of G. The restriction K~ i ^ c of K 7 to the reduced 
phase space corresponding to the value c of G is then given by 

K^ C {L X ,4>) = i (d!, 7 (r 2 - I 2 ) + d 2n (r 2 2 - {Lx - c) 2 )) 

+ ^{r 2 -L 2 1 ){r 2 -(L 1 -c) 2 ) cos0 (65) 

with L x S ((-n,n) n (c- r 2 ,c + r 2 )) and G M/2ttZ. 

The reduced Hamiltonian vector field induced by -ftT 7 , c is given by 

d f L A _ fr „u f 9K-r,c/0<l> 

Note that 

{ il ,0 } = {i 1 ,axcta n ^} = T - R L-^{i 1 ,l} 
1 ct{£i,t} - rji^cr} 



x^z^) = dt{t)= {L ^ } { -bkZ/bl, ) • (66) 



1 + (t/ct) 2 a 2 l + (r/(i) 2 cr 2 

since {Xi,t} = cr and {Li,er} = — r. Furthermore, with L 2 = L\ — c, 



^ - -'(^-^(rl-Ll) sin, 



-— — - ~{d l , 1 L 1 +d 2 ,- l L 2 ) _______ C os < 

Hence (1661) reads 



^(ir^) = (d Li + d L2) + ^f^^) C os J' ( 67 ) 
where we treat L 2 = £1 — c as a dependent variable. 
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5 CHAINS WITH DIRICHLET BOUNDARY CONDITIONS 



By (1571) . the fixed points of the vector field X JtC with \Li\ < r\ are given by 
(L\, 4>) satisfying 

(j) G ttZ (68) 

and 

(4, 7 £! + d 2 , 7 i 2 ) + ^fa^^fe^ cos * = 0. (69) 

Note that |68|) and the square of ([69]) are equivalent to the system of equations 
(j52" j) -([59 ]) derived above. 

In order to determine the type of the fixed points (Li,</>), i.e. points sat- 
isfying (|68|) - (|69|) for a given value c of G, one computes the Jacobian of X 1}C , 

H = Hry C (Li,<b) — ( (* 12 ] at these points. Note that at such points 

V "21 "22 / 

tin = and /122 = and thus det(H) = —hx^h^i- Hence such a fixed point 
is an elliptic or hyperbolic fixed point of X ltC if /112/121 is positive, respectively 
negative. We omit a more detailed analysis of these points. 



5 Chains with Dirichlet boundary conditions 

In this section we consider a FPU chain with N' (N' > 3, not necessarily even) 
moving particles and fixed endpoints, i.e. with boundary conditions (|18[) . 

It has been observed that such a chain can be treated as an invariant sub- 
sytem of a periodic lattice with N = 2N' + 2 particles - see [13]: Let T*R N be 
endowed with the canonical symplectic structure and consider the linear map 
S : T*R N — > T*R N mapping (qi)i<i<N, (Pi)i<i<N to 

-(qw-i, ■ ■ -,qi,qN), -{pn-1, ■ ■ ■ ,Pl,PN)- 

Then S is a canonical linear involution satisfying Hy o S = Hy. Denote by 
Fix(S') the fixed point set of S. Then Fix(S') is the subset of all elements (q,p) 
in T*R N satisfying 

(q n ,Pn) = -(qN-n,PN-n) VI<n<iV-I and q N = p N = 0. (70) 

In particular, on Fix(5), qN — qN'+i — and pn = Pn'+i = 0. Note that on 
Fix(S'), both the center of mass coordinate Q = -h 53-t=i ft an< ^ its momentum 
P = J2iLiPi are identically 0. Hence Fix(S) C M, where 

M := {(q,p) G T*R N \Q = 0; P = 0}. 

We endow M with the symplectic structure induced from T*R N . 

The phase space of an FPU chain with N' moving particles satisfying Dirich- 
let boundary conditions is T*R N , endowed with the canonical symplectic struc- 
ture J2iLi dqi A dpi. It can be embedded into M by the map 6 : T*R N ' — > M 
defined by 



(ft,Pi)l<KJV' !-* —7= {(qi,Pi)l<i<N', (0, 0), -(qN'^i,PN'^i)o<i<N'-l, (0,0)) 
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Note that 9(T*R W ) = Fix(S), i.e. 6 is a parametrization of Fix(S) and 

the pullback of the canonical symplectic form on M by 9 is J2i=i d 1i A d Pi> 
which means that 9 is canonical. It then follows that Fix(S') is a symplectic 
submanifold of M. 

We now express the equations defining Fix(S) locally near as a subset of 
M in terms of the canonical coordinates (xk, Uk)i<k<N-i provided by Theorem 
II. 1[ or even more conveniently, in terms of the associated complex coordinates 
(Ck)i<|A|<J\r-i) defined for 1 < k < N - 1 by 

J a = ^(a*-Hfc) ( 

\ C-fc = C* = -^(xk+iyk). [ } 

Denote by Z the linear subspace of C 2N ~ 2 consisting of such vectors ( Cfe ) i<| fc |<iv— i - 
In the sequel we also write (Ck)i<k<N-i for the element (Cfc)i<|fc|<jv-i G Z and 
use the notations (n £ Z) 

n7T T17T 

c n := cos — , s„ := sin — . 
Define the map Sz '■ Z —* Z, given by 

(Cfc)l<ft<JV-l (-e imk l N (N-k)l<k<N-l- (72) 

Like the map S : M — > M, Sz is a canonical linear involution. In fact, the 
maps S and Sz are conjugate to each other under the coordinate change of 
Theorem 11.11 Before making this statement more precise, let us introduce a 
parametrization of the fixed point set Fix(Sz) of the map Sz- Introduce 

Zn ir ■= {(00i<|*|<jv e C 2N '\( k = C-k VI < k < N'}, 

endowed with the canonical symplectic structure induced from C 2N , and the 
embedding Q z ■ Z Dir — > Z mapping (Cfc)i<|fc|<JV to the element (Cfc)i<fc<jv € Z 
given by 

— ^ ^(Cfc)i<fe<jv, 0, (— e 4,rlfc / A 'CA r '+i-fc)i<fc<Af'^ • 

Note that Oz(Zz?ir) = Fix(Sz), i.e. Gz is a parametrization of Fix(Sz). 

Using the explicit construction in [7] of the canonical transformation ^ of 
Theorem 11.11 one verifies the following lemma. 

Lemma 5.1. In terms of the complex variables (Cfe)i< |fe|<JV l defined by The- 

orem \Tj[ nearO, the map S is given by Sz- More precisely, if ^ , defined near 
G Z, is the coordinate transformation given by Theorem then S o = 
o Sz- In particular, locally near 0, the set Fix(Sz) C Z, described by the 
equations 

e- 2mk/N ( k + e 2mk/N QN-k = (1 < k < N - 1), (73) 
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5 CHAINS WITH DIRICHLET BOUNDARY CONDITIONS 



is the image of the set Fix(S) under '5 . Expressed in terms of the real variables 
(xk,yk)i<k<N-i, the conditions (731) are given by 



C2k ~S2k \ f Xk \ , f C 2 k S 2 k \ ( XN-k \ _ f 
S2k C 2 k J \Uk J \ ~S 2 k C 2 k J \ VN-k ) \ 



(74) 



In particular, for k = N' + 1 (= N/2) we get CiV'+i — and therefore 

(xn'+i,Vn'+i) = (0,0). 
Corollary 5.2. On Fix{S z ), for anyl<k<^, 

h = I N -k (75) 

and 

J k J § _ k - M k M § _ k = hlN_ k . (76) 

Moreover 

If = 0. (77) 

Proof. In terms of the complex variables (Cfe)i<|fe|<iv-l; Ik = CfcC-fc f° r an y 
1 < k < N- 1. Hence on Fix(Sz), 

h = CfcC-fc = (-e 47rife / Ar Cjv-fc)(-e _47rlfc ^ A 'C-(Af-fc)) = CiV-fcC-(JV-fc) = 

showing ([TBI . The identity (f77|) follows from C^|Fte(S z ) = 0. To prove (f76|) . we 
first conclude from ((73)) that on Fix(S*z), for any 1 < k < iV — 1, 

«/fe = —C4klk, 
M k = -s 4 fc7fc. 

Hence on Fix(Sz), 

JkJf^u - M k MiL_ k = /fc/«_ fe (c4fec 4( «_ fc) - s 4fc s 4( «_ fe) ) 

= 4-T*_fc(C4fc + s 4fc) 
= Ikl^-h- 

2 

This completes the proof of Corollarv l5.21 □ 

From the definitions ([6]), (jlip . and (fl2|l of the variables Jfc, Mfc, and of 
the expressions R and Rn_ one then obtains the following 



Corollary 5.3. On Fix{S z ), 

R = 4 S2klkl <Mld i?/ 



l<fe<^ 



4 «/fe 

4 

otherwi 



It follows from Corollary 15.31 that on Fix(5z), the expression (fT3)) is in 
Birkhoff normal form up to order 4. This allows us to prove Theorem 11.31 
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Proof of Theorem \1.3l We start with the resonant normal form p3[) for even 
chains, ^-+H a ^(I)~R a ^(J, M)+0(\(x, y)\ 5 ), where H a ,p(I) and R a ,p{J, M) 
are given by ([5J) and (jTDJ) , respectively. Using the identity If. = In-Hi the terms 
in the decomposition Q14p of H a p, when restricted to Fix(5 l 2,), are given by 



N' 



fc=i 



N' 



fe=l 



4(1 -a 2 ) 
27V 



^ SkSihk- 



and 



l<fc,I<itf' 



AT' 



1 % 1 \ 



2N 



(78) 
(79) 

(80) 



k=l k=l 

From Corollary 15.31 we conclude that on Fix(5z), 

[3 -a 2 



AN 



(3 -a 2 
4iV 



(r(J,M) +Rn(J,M)^ 



( 



\ 



4 ]T S2khlN_ k + 1% 

only if f-ew/ 



Kfc<^ 



(81) 



Formula (fTi?)) is then obtained by adding up (|78 p -(|8i p . noting that dt + -£- = 



\{a 2 + 3(/3 — a 2 )s^), and replacing Ik by its pullback with respect to the 
parametrization Qz of Fhc(Sz) introduced above. □ 

It remains to prove Theorem 11.51 We first consider the case a — 0. 

Proposition 5.4. Assume that a — in Then the following holds: 

(i) The Birkhoff normal form of Hy with Dirichlet boundary conditions up to 



{N' + 1)P 2 



order 4 is given by > 

N' 



+ Hgp(T) where 

N' 



N ft f 1 

k=l y ' \k=l 



only if ^tiei 



N> \ 

"4 ^ SlS m IlI m — ^ S2kTklN' + l-k ) 
M"> fc = i / 



(82) 



l<i,m< JV 7 



ftij For any (3^0, H^g(I) is nondegenerate at I = 0. 
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5 CHAINS WITH DIRICHLET BOUNDARY CONDITIONS 



Proof. The Birkhoff normal form (|82[) of Hy with Dirichlct boundary conditions 
is given by the formula (fTO)) evaluated at a = 0. To investigate the Hessian of 



Q° of H^g(J) at I = 0, we write 



16(JV' + 1) 



A N'pD A N> 



(83) 



where A w = diag ( sin „ „ ) and P 15 is the N' x iV'-matrix which for 

V ZN +z /l<k<N' 

N' even resp. odd is of the form 



/3 4 ... 

4 3 4 



... 4 2\ 

4 2 4 



3 4 4 2 

4 3 2 4 
4 2 3 4 
2 4 4 3 



4 2 4 
\2 4 ... 



4 3 4 
...4 3/ 



/3 4 ... 

4 3 4 



... 4 2\ 

4 2 4 



3 4 2 

4 2 4 
2 4 3 



4 2 4 
\2 4 ... 



4 3 4 
...4 3/ 



(N' even) 



(N' odd) 



where we used that S2k = 2stCk = 2skSN>+i-k and, if jV 2 +1 e N, s^ /+1 

2 

follows that 



It 



det Q 



D 

0,0 



2/3 



16(iV' + l) 



AT-1 



N' 



• det P° • Y[ sin 2 



fc/T 



fc=l 



2iV' 



In order to see that P D is nonsingular, observe that detP D S Z. For A/ 7 even 
we show that detP-° = 1 mod 2. Note that in this case the diagonal of P D 
consists of 3's only. Therefore det P D = "S N mod 2 = 1 mod 2. If N' is odd, the 
same argument shows that det P = 2 mod 4. Hence, if [3 ^ 0, det Qj^g ^ 0, and 
the nondegeneracy of the Hessian of HqJI) at I = follows. □ 

Lemma 5.5. If (3 < 0, then Qqo has r JV 2 +1 ~ I negative eigenvalues, whereas if 



N'-l 



(3 > 0, then Qq p has l ,, 



negative eigenvalues. In particular, for any (3 =/= 0, 



q ^ is indefinite (and HP ^ is therefore not convex). 



Proof. We want to use the decomposition 



be deformed continuously to 



2/3 
16(W'+1) 



(*A iV ' + (l-t)Id) (tA iV ' +(!-*) Id) 



of Qj^ to show that Qq can 



P D : Consider for < t < 1 



>AT' 



16(JV' + 1) 
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As t A w + (1 — t) Id is positive definite for any < t < 1 and P D is regular and 
symmetric, Qq «(t) is a symmetric regular N' x TV'-matrix for any < t < 1. 
For t = 0, Q^'(0) = 16( ^ +1) P D , whereas for * = 1, Q^ (3 (l) = Q^. Therefore, 
index(Q^ > /3 ) (i.e. the number of negative eigenvalues of Qq a) coincides with 
index( 16 r^f +1 s P D )- To describe the spectrum of P D , we distinguish between 
N' even and odd. 

If N' is even, the eigenvalues of P D are 47V' — 3 (with multiplicity one), 1 
(with multiplicity ^-), and —3 (with multiplicity ^ — 1), hence P D has ^ — 1 
negative eigenvalues. If N' is odd, the eigenvalues of P D are 1 (with multiplicity 

^=i), -3 (with multiplicity £^5), and i(4/V'-5) (l ± yjl + ff^j l) (each 

with multiplicity one), hence P D has jV 2 ~ 1 negative eigenvalues. These facts 
are verified in Appendix [Al The claim of the lemma now follows. □ 

We now turn to the case a^O. 

Proposition 5.6. Assume that a ^ in Then, for a fixed, detQ^ p * s a 
polynomial in (3 of degree N' and has N' real zeroes (counted with multiplicities) . 
When denoted by flk = flk{ot) (1 < k < N' ) and listed in increasing order, they 
satisfy 

Pi< ■■■< PriH+in < a 2 < l3 rEl +3^ <■■■< Pn>- 

2 2 

Moreover 

for [3<f3i 
for [3 

2 2 

/or /3 > /3jv 

Proof. Fix a e K \ {0} and consider the map (3 i— > det(Q^ ) ^). It follows from 
(fl9|) that det(Q^ ) ;3 ) is a polynomial in /3 of degree at most N', 

N' 
3=0 

where r = det(Qf ) and rjv = det(Q^ 1 ). By Proposition [Ql det(Q^ 3 1 ) ^ 0, 
hence the degree of the polynomial det(Q® ^) is TV'. We claim that det(Q^) has 
N' real zeroes (counted with multiplicities). For |/3| large enough, index(Q^ ) /3 ) 
is equal to index(Q^). By Lemma f5.5[ index(Q^) is l-^— -j for [3 > and 
r jv^fin for ^ < o. Hence there exists R > such that index(Q^ /3 ) = L^^ij 
for any /3 > R and mdex(Q£ i/3 ) = ri ^ tin for any f3 < -R. For (3 = a 2 , 
Qa a 2 ^ S a positive multiple of the identity matrix, hence index(Q^ Q2 ) = 0. It 
then follows that, when counted with multiplicities, index(Q^) must change 

at least r jv 2 +1 ~ l times in the open interval (— oo, a 2 ) and at least i_ jv 2 ~ 1 _i times 
in (a 2 , oo). Since a change of index(Q^ ) /3 ) induces a real zero of det(Q^ ) (g ), our 



index (Q^p) 



r jV_±l -| 

2 
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A SPECTRUM OF THE MATRIX P 



consideration shows that (3 i— ► det(Q® (3 ) has N' real zeroes. Further we have 
/3 r jv' + i -i(g) < a 2 < /3 r iy/ + 3 -|(a). This proves the proposition. □ 

2 2 

Proof of Theorem \1.5[ Part (i) is proved by Proposition l5.6l whereas (ii) follows 
from Propostion l5.4l and Lemma T5. 51 □ 



A Spectrum of the matrix P 



D 



Here we compute for any integer N' > 3 the eigenvalues of the N' x N' matrix 
P D , given by 



/3 4 ... 

4 3 4 



4 2 4 
\2 4 ... 



... 4 2\ 

4 2 4 



3 4 4 2 

4 3 2 4 
4 2 3 4 
2 4 4 3 



4 3 4 
...4 3/ 



(N' even) 



(2, 4 ... 

4 3 4 



3 4 2 

4 2 4 
2 4 3 



4 2 4 
\2 4 ... 



... 4 2\ 

4 2 4 



4 3 4 
...4 3/ 



(jv odd) 



Lemma A.l. If N' is even, the eigenvalues of P D are AN' — 3 (with multiplicity 

one), 1 (with multiplicity ^-), and —3 (with multiplicity ^ lj. J/ TV' is odd, 

the eigenvalues of P D are 1 (with multiplicity ^ )> ~3 (with multiplicity 
jV 2 ~ 3 ,), and i(4Af' — 5) f 1 ± ^/l + || ^P~^p ^ (each with multiplicity one). 

Proof. Throughout this proof, antidiag(ai, . . . , ajv') denotes the "antidiagonal" 
AT' x TV'-matrix M with Mjy = a ; if fc + I = N' + 1 and M fe; = otherwise. 
First consider the case where N' is even. We write P D in the form 

P D = diag(-l, . . . , -1) + antidiag(-2, . . . , -2) + 4 • l N , xN , 

and, with /i := —1 — A, 

P D - Aid = diag(/i, . . . , /z) + antidiag(-2, . . . , -2) +4 ■ 1 



N'xN'- 



= :L(«') 



Here Iat'xw denotes the N' x Af '-matrix whose entries are all equal to 1. 

We compute dct(P D - Aid) = detfX^'' +4-1 N'xN') by column expan- 
sion. Note that in the column expansion of the determinant only those terms 



31 



contribute which are determinants of matrices containing at most one column 
consisting of entries all equal to four. We obtain 

N' 

det(P D - Aid) = det(L^ N ">) + J] det^f' 3 ), (84) 

3=1 

where Lj is defined as the matrix L^ N ■* with the j-th column replaced by the 
column 4 • ljvxi- By expansion with respect to the first column and then the 
last column, the determinant of L^ N ' can be computed recursively, 

det(i( JV ')) = ( / i 2 -2 2 )det(i( jv '- 2 )) 

Since det(L^ 2 )) = /i 2 — 4, it follows by induction that 

det(L( JV ')) = (/x 2 -4)T. (85) 

To compute det(L[ N '), we expand the determinant in the same way and obtain 
the identity de^L^) = 4(/x + 2) det^'- 2 )), fr om which it follows that 

det(L[ N ">) = 4(/x + 2)0? - (86) 

Similarly one gets det(L^ "*) = (p 2 — 4) det(L^ 2 '), and thus 

det(4 JV ' ) ) = det(L{ N ' } ) = 4( M + 2)(/i 2 - 4)^ i - 1 . (87) 

For any 1 < j < this procedure leads to det(ij JV ') = {p 2 — 4) det(L^ 1 -2 ^) 
and hence 

det^ ) = det{L[ N ' ] ) = 40* + 2)(^ 2 - 4)*- 1 . (88) 

For j > 4r, note that detL^- 1 = detL^^-, since Lj N ^ and L^Jj can be 
transformed into each other by exchanging the j'th and the (N — j + l)'th 
columns and then the j'th and the (N — j + l)'th rows. By (|85 p -(|88 p . we obtain 

det(P D - Aid) = (fi 2 - 4)^- x • (0* 2 -4) + N'- Mp + 2)) 
= On 2 - 4)^- 1 ( A t + 2)(/i - 2 + 47V'). 

Hence, if N' is even, the zeroes of det(P D — Aid) are /i = 2 (with multiplicity 
T! — 1), /i = —2 (with multiplicity ^-), and /i = — 4iV' + 2 (with multiplicity 
1). Transforming back to A = —1 — /i, we obtain the claimed eigenvalues. 
It remains to consider the case where N' is odd. Again, we write 

P D =diag(-l,...,-l, ^ ,-l,...,-l) + antidiag(-2,...,-2)+4-lAr/ xA r/. 
With [a = — 1 — A we get 

P D - Xld = L^ N "> + 4- l N/xN , 
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with 



/ H o 

o '■• 



L (N> 



-2\ 

■ ■' 



1 



o ••' '■• 

V-2 ... \i ) 

As above, we obtain the expansion ([84]) for the determinant of P D — Aid. We 
expand dct(L^ Ar ') with respect to the first column and then the last column, 
yielding the recursion formula 

det( J L( JV ')) = (/i 2 -4) det(L( w '- 2 >), 

which together with det(L^ 1 )) = jj, — 1 leads to 



det{L^) = {^- 4)^-1). 



(89) 



For det(L[ N ">), we obtain the identity det{L[ N ' } ) = 4(/i + 2) det(L^ N '-^) and 



hence 



det(L[ N,) ) = 4( M + 2)( M 2 - 4)^ (/i - 1). 

N' 



More generally, for any 1 < j < 4-, we have 



detiE^detCif'^det^). 



(90) 



(91) 



It remains to compute det L . Expanding det L by the first column and 

2 2 

then the last column, we obtain the recursion relation 

det L^'l = (fj 2 - 4) det L { ( ^ +1 . 



Together with det L 2 




4(/i 2 — 4), this implies 



(92) 



deti£2 i =4( M 2 -4)^. 

2 

Hence, combining (f89|) - (|92|) we obtain 

det{P D - Aid) = (a^ 2 -4)^ • (( M 2 -4)(/^-l) + (iV'-l) • 4(/i+2)(/i-l) + 4( M 2 -4)) 

= (^ 2 -4)^0*+ 2)(^ 2 + {AN' - 3)m - {AN' + 2)). 
Hence, if N is odd, the zeroes of det{P D — Aid) are ji = 2 (with multiplicity 



N'-3 



), /i = — 2 (with multiplicity 



N'-l 



), and 



// = --{AN' - 3) ± i\/l6iV' 2 -8^' + 17 

(each with multiplicity 1). Transforming back to A = —1 — /i, wc obtain the 
claimed formulas for the eigenvalues in the case where N' is odd. This completes 
the proof of Lemma IA.ll □ 
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